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Abstract 

We discuss the Hamiltonian H = — (ix)^"^^ and the mixed Hamiltonian 



^ ■ -f^mixed = (p^+x^)/2 — (^(ix)^"^^ . The Hamiltoniaus H and in some cases also -ffmixed 

are crypto-Hermitian in a sense that, in spite of their apparent non-Hermiticity, a 
quantum spectral problem can be formulated such that the spectrum is real. We 
0^' note that the corresponding classical Hamiltonian system can be treated as a gauge 

system, with imaginary part of the Hamiltonian playing the role of the first class 
^ i constraint. 

^ ' Several different nontrivial quantum problems can be formulated on the basis of 

this classical problem. We formulate and solve some such problems. We consider 
then the mixed Hamiltonian and find that its spectrum undergoes in certain cases 
a rather amazing transformation when the coupling g is sent to zero. There is an 



(7) 

' infinite set of exceptional points gl where a couple of eigenstates of H coalesce and 

. their eigenvalues cease to be real. 



When quantization is done in the most natural way such that gauge constraints 
are imposed on quantum states, the spectrum should not be positive definite, but 



Q ■ must involve the negative energy states ( ghosts). We speculate that, in spite of the 

, appearance of ghost states, unitarity might still be preserved. 

O 



1 Introduction 

For certain apparently complex Hamiltonians, the spectral problem can be formulated 
such that the spectrum has a perfectly "normal" form with bounded from below real 
energies. Such Hamiltonians can thus be called "crypto-Hermitian" or "cryptoreal" . 
Apparently, such crypto-Hermitian Hamiltonians were first discussed in association with 
Reggeon field theory back in 1976 [1]. Somewhat later, crypto-Hermitian Hamiltonians 
were considered by mathematicians in a more habitual Schrodinger setup. Gasymov 
observed that the Schrodinger operator with certain complex periodic potentials, like 
V{x) = e*^, has a real spectrum [2]. In Ref. [3], it was proved that the spectrum of the 
Hamiltonian with complex potential V{x) = x'^ + i(3x^ is real for small enough (3. General 
properties of crypto-Hermitian (or quasi- Hermitian as the authors called this property) 
operators were studied in Ref. [4]. 



*0n leave of absence from ITEP, Moscow, Russia. 
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Before going further, a comment on the terminology is in order. At the moment, 
there is no unique generally adopted name for Hamiltonians of this kind. Besides quasi- 
Hermitian, the term pseudo-Hermitian is also often used. A Hamiltonian is usually called 
pseudo-Hermitian if it satisfies the property 



with some Hermitian invertible rj. However, this does not guarantee yet that the spectrum 
is real. To this end, the operator rj should be representable as [5] 



or, equivalently, the norm (-0 \r)ip) should be positive definite for any nonzero Hilbert 
space vector ip [6], [4]. Anyway, the semantics of the words quasi-Hermitian or pseudo- 
Hermitian is "not quite Hermitian" with a flavour of inferiority, "second-rankness" com- 
pared to Hermitian. For example, pions are pseudo-Goldstone particles meaning that 
they are not Goldstone particles. But we want to emphasize that, if the spectrum of the 
Hamiltonian is real, the latter almost always is in fact Hermitian when looking at it 
through proper glasses, i.e. when defining the norm in Hilbert space in a proper way. It 
was proved in [5] that the Hamiltonian with real non-degenerate spectrum must satisfy 
the properties ([H [2]). Then r] defines the norm with respect to which the Hamiltonian H 
is Hermitian, while the Hamiltonian H = 0H0~^ is manifestly Hermitian with respect to 
the standard norm. In other words, the characterization "crypto-Hermitian" (Hermitian 
in disguise) reflects more adequately, in our opinion, the essence of the phenomenon, and 
we will stick to it in this paper. [f| 

The modern history begins with the beautiful paper [9] (see also the recent review [10]), 
where this property was observed for a wide class of PT-symmetric polynomial potentials, 
like V{x) = ix^. 0. It was found to be discrete and real. 

Since then, many crypto-Hermitian Hamiltonians have been discovered. We can men- 
tion the paper [11] where the spectrum of the Hamiltonian with hyperbolic and generalized 
hyperbolic complex PT-symmetric potentials was shown to be real in many cases. The 
simplest example of such cryptoreal hyperbolic problem is the problem with the potential 



with Vi > and IV2I < Vi + 1/4. 

In recent [12], it was shown that apparently complex Hamiltonians obtained after so 
called nonanticommutative deformations [13] of certain supersymmetric quantum-mechanical 
and field theory models are in fact crypto-Hermitian and enjoy a real spectrum. 

^ "Almost" means away from exceptional points [7] where the Hamiltonian involves Jordan blocks. We 
will discuss this issue later. 
^It was used first in Ref. [8] 

^Let us repeat for clarity: our crypto-Hermiticity means exactly the same as quasi- Hermiticity of 
Ref. [4] (a quasi-Hermitian Hamiltonian was defined there as the Hamiltonian that is Hermitian with 
respect to a generalized positive definite norm {ip\riip), rj^ — rj), the same as rj- Hermiticity as defined in 
Ref. [6] and the same as pseudo- Hermiticity ([T]) with additional requirement ([2|). 
A PT-symmetric potential V{x) enjoys the property V*{—x) = V{x). 



H^ = rfHrj 



-1 



(1) 



r] = 0^0 



(2) 




(3) 
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The problems with the potential V{x) = e*^ or the potential ([3]) admit explicit analytic 
solutions. In [9], reality of the spectrum for the potentials V{x) = x'^iix)" , e > 0, was 
demonstrated explicitly by numerical solution of the corresponding Schrodinger equations 
supplemented by semiclassical analysis. Later, a rigorous proof for the discreteness and 
reality of the spectrum in this problem was constructed [14]. In Ref. [15], it was shown 
that the Hamiltonians like 

H = - g{^xr-^' (4) 

can be represented for small g in the form ([ll [2]). In other words, they can be obtained by 
a non-unitary transformation, H = e~^He^, out of a manifestly Hermitian Hamiltonian 
H. The Hamiltonian H and the operator R = InO are calculated perturbatively as an 
infinite series in the coupling constant g. 

In this paper we suggest an approach capitalizing on a certain hidden gauge sym- 
metry characteristic of crypto-Hermitian systems. The origin of this symmetry is very 
simple [16]- [18]. Consider a system with one dynamical degree of freedom. The classical 
Hamiltonian is a function x), which may be real or complex. Let us complexify the 
phase space variables, 

X ^ z = X + iy, p ^ n = p — iq , 

n{p,x) ^n{n,z) = H{p,q;x,y) + iG{p,q;x,y) , (5) 

where H and G are real functions satisfying the Cauchy-Riemann relations 

dH dG _dH dG _ dH dG _dH dG _ ^ 
dp dq dy dx dq dp dx dy 

Two important properties follow : 

• The function H{p^ q; x, y) can be treated as the Hamiltonian of a new system with 
double set of degrees of freedom. Indeed, the real and imaginary parts of the com- 
plexified equations of motion for the original system, 

TT = —dH/dz, z = dH/dn , 

coincide in virtue of ([6]) with the Hamilton equations of motion derived from H{p, q; x, y) . 

• The Poisson bracket 

{H G} , - + dl^OG _ dl^dG _ dH dG 

dx dp dy dq dp dx dq dy 

vanishes. This means that G is an integral of motion for the system described by 
H. The space of all classical trajectories is thus divided into classes characterized 
by a definite value of G. The class with G = represents a particular interest. The 
condition G = can be interpreted as a first class constraint and the dynamical 
system with the Hamiltonian H supplemented by the constraint G = is a gauge 
system. 
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The plan of the paper is the following. In the next section, we consider from this angle 
the simplest possible problem — the complexified oscillator. We note that this classical 
problem has at least three different quantum counterparts: 

1. One can impose the analyticity constraint on the wave function, dip/dz = and 
solve the Schrodinger equation in the vicinity of the real axis. In this case, we 
reproduce the standard oscillator spectrum En = 1/2 + n. The same spectral 
problem is obtained when the gauge constraint is resolved at the classical level with 
the gauge choice y = 0. 

2. One can impose the analyticity constraint and solve the Schrodinger equation in the 
vicinity of the imaginary axis. In this case, the spectrum En = —1/2 — n involves 
negative energies and is bounded from above rather from below. The same spectral 
problem is obtained when the gauge constraint is resolved at the classical level with 
the gauge choice x = 0. 

3. Finally, one may not require analyticity, but rather impose, following Dirac, the 
gauge constraint G"^ = on quantum states. In this case, the spectrum is En = n, 
where n can be positive, zero, or negative. Still, the quantum problem is well 
defined, and the evolution operator is unitary. 

In Sect. 3, we consider the classical dynamics of the Hamiltonian 



We find different sets of trajectories with positive, and also with negative energies. Sect. 
4 is devoted to the quantum dynamics of ([7]) and of the mixed Hamiltonian (jl]), with 
the analyticity constraint imposed on wave functions. The complex plane of z is divided 
then into several regions. In some of them, the spectrum is discrete, in some others - 
continuous or empty. For n = 1, we reproduce the results of Ref. [9]. For n > 1, one can 
formulate n different spectral problems with discrete positive definite spectrum formulated 
in the different regions of the complex plane. When g ^ 0, the spectrum of the mixed 
Hamiltonian approaches the oscillator spectrum, but for the problems formulated in the 
sectors not comprising real axis, the transformation pattern is very nontrivial involving 
an infinite set of "phase transitions" in the coupling. At each of such "phase transition" 
(or exceptional [7] point gi''^), a pair of eigenstates of the mixed Hamiltonian coalesce 
such that at this very point the Hamiltonian involves a nondiagonalizable Jordan block. 
Beyond this point {g < gi''^), a pair of complex conjugate eigenvalues should appear. In 
other words, the Hamiltonian is cryptoreal in this situation only for large enough 
g (and, obviously, it is Hermitian for g = 0). The last section is devoted as usual to 
discussions. In particular, we discuss the Dirac spectral problem for the Hamiltonian ([7j) 
when the gauge constraint is imposed on the wave functions as an operator condition. 
This problem has no analytic solution and is difficult to resolve numerically. Still, based 

^That was observed in Ref. [19]. The proof of reality and discreteness of the spectrum for aU symmetric 
(see below) spectral problems was given in Ref. [20]. 




(7) 
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on semiclassical reasoning, we argue that, similar to what we had in the oscillator case, 
the spectrum there might be discrete and unbounded both from below and above. We also 
point out the similarity of this problem to some other previously analyzed by us systems, 
which are described by higher derivative Lagrangians and involve ghosts. We speculate 
that, in spite of their presence, unitarity is not violated. 

2 Complex oscillator. 

Consider the complex Hamiltonian 

7i(7r,^)=^^. (8) 

Its real and imaginary parts are 

H = l^-l±>^ (9) 

and 

G — —pq + xy . (10) 
Consider the classical dynamics of H. The classical trajectories are 

X = Asm.it + p = Acosit + 4>i), y = B sm{t + 4>2) , q = —B cos{t + 4>2) ■ (H) 

Generically, they have complex energies. If we require the energies to be real, i.e. impose 
the constraint G — 0, the relation 

cos((/)i - 02) = (12) 

follows. For each value of the energy, positive or negative, there is a set of trajectories 
(cofocal eUipses) with the same period (see Fig. 1). 

In the case under consideration, the period is the same for all energies, but this is the 
specifics of oscillator. The fact that the period is the same for all trajectories of a given 
energy has, however, a general nature. In fact, it is a consequence of the gauge symmetry 
of the problem. 

The latter is simply the symmetry generated by the constraint G. Infinitesimally [18], 

5gx = -a{G, x}p,B. = -aq, Say = -a{G, y}p.B. - ap, 

5gP = -a{G,p}p,B. ^ -ay, Saq ^ -a{G,q}p.B. ^ -ax . (13) 

This is a phase space symmetry. To represent it as a conventional gauge symmetry acting 
only on the coordinates, one should introduce the Lagrange multiplier X{t) and write the 
canonical Lagrangian as 

L = px + qy-H-XG , (14) 
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Figure 1: Family of oscillator trajectories with the same positive energy. For negative 
energies, the picture is rotated by 7r/2. 



Expressing out the momenta, 



x-Xij y + Xx 



we obtain 



^= 2(i + A^) +^^->^^y- (16) 

The gauge transformations amount to shifting the Lagrange multiplier A by (a derivative 
of) an arbitrary function of time a{t), supplemented by the transformations of dynamic 
variables x, y generated by the constraint G. 



a{y + Xx) 

gauge-^ — — Lxy — ^ _^ ^2 ' 

a{x — Xy) 



^ gauge O^Q 

Sgaugey = —<yp 1 + ' 

^gaugeA = a . (17) 

Indeed, one can explicitly verify that the Lagrangian f|T6l) is invariant, up to a total 
derivative, with respect to the transformations f|T7|) . 

The transformations Sx and 6y in Eqs. ( fT3l [T7l) have a clear meaning. Any Hamiltonian 
system is invariant with respect to time translations t ^ t — a that transform a solution 
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z{t) z(t — a). Their generator is the Hamiltonian H. In our case, however, besides 
H = Re(7^), we have another integral of motion G = Im(7^). It generates a shift of time 
by an imaginary amount, t ^t — ia and transforms z{t) z{t — ia). Infinitesimally, this 
coincides with Eq. flTTI) (with partial gauge fixing A = 0). 

The shift z{t) z(t — a) is the shift along the trajectory, leaving it unchanged. But 
the shift z{t) — > z{t — ia) transforms one trajectory into another. It is this shift which 
relates different ellipses in Fig. 1 [it is straightforward to check by substituting for z{t) 
the exact analytic solution ffTTj) with 01 — 02 = 7r/2]. Such families of closed trajectories of 
a given energy and the same period (obviously, if z(t) is periodic, z(t — ia) is also periodic 
with the same real period) exist also for more complicated cases. We will discuss it in the 
next section. 

Let us go over to quantum dynamics. There are two basic ways to quantize gauge 
systems!^: (i)hy explicitly resolving the constraints and quantizing the Hamiltonian with 
a reduced number of degrees of freedom; (ii) by not resolving the constraints classically, 
but rather solving the system 

H^f = E^, = . (18) 

We will see that, in the case under consideration, these two approaches are not quite 
equivalent, in contrast to what is usually assumed ! 

• Let us first try to resolve the constraint G = classically. This can be done by 
fixing the gauge, i.e. by imposing the additional constraint xip-il'^^TU) = 0, where 
{G, x}p.B. 7^ (so that the primary constraint G = and the gauge fixing constraint 
X = are independent). Resolving the system G = x = 0, we are left with a reduced 
number of dynamical variables. Generically, their number is equal to the number 
of initial degrees of freedom minus the number of primary constraints. In our case, 
-/Vreduced = 2 — 1 = L One can, for example, choose x = 1/ = 0. The reduced 
Hamiltonian system will in this case be just H* = {p^ + x^)/2 with the spectrum 
En = 1/2 + n. On the other hand, if choosing the gauge x = x = 0, the reduced 
Hamiltonian is H* = —{q^ + l/^)/2 with a different spectrum = —1/2 — n. In 
other words, there are two essentially different gauge choices leading to different 
reduced Hamiltonians. One can obtain either oscillator with positive energies, or 
oscillator with negative energies, but not both. 

To understand what happened, look again at the trajectories in Fig. 1. They 
represent, as we have seen, gauge copies of one another. The gauge fixing procedure 
should pick out one of these copies, while getting rid of all others. And, indeed, 
the condition y = does this job by pinpointing the trajectory going along the real 
axis. However, none of these trajectories are compatible with the condition a; = 0. 
On the other hand, for the family of the trajectories with negative energies, one can 
impose a; = (and pinpoint the trajectory going along the imaginary axis), but not 

^This problem was first posed and resolved by Dirac and is treated pedagogically in many books. See 
e.g. Ref. [21] 
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?/ = o.0 

Two spectral problems with positive and negative energies can alternatively be 
defined using the approach of Ref. [9]. To this end, one should require that the wave 
function represents an analytic function of z = x + iy. The spectrum En = 1/2 + n is 
then realized by the standard oscillator functions continued analytically to complex 
arguments. For example, the wave function of the ground state is exp(— 
It falls down exponentially on the real axis and also on the lines z = us, s G 
(—00, 00), |Arg(M)| < 7r/4. The spectrum En = — 1/2— n is realized by the functions 
like exp(z^/2) that fall down exponentially along the imaginary axis and in the sector 
|Arg(z)| > 71 /A (see Fig. 2). 











E<0 































Figure 2: Sectors with positive and negative energies in the spectrum 

The lines Arg(z) = ±7r/4 are closely related with the Stokes lines of the oscillator 
Schrodinger equation. The Stokes lines are defined [22] as the lines that pass through 
turning points and satisfy the condition 

Im (^j^ Tr{w)d'u?j = (^j ^ ^2{E - w'^)dn^ = , (19) 

( zo is the position of the turning point). The asymptotes of Stokes lines at large 
values of \z\ are the straight lines separating the sectors in Fig. 2. When crossing a 



""The trajectories in Fig. 1 are related by gauge transformations with constant a. But one can easily 
prove that one cannot obtain a configuration with x(t) = out of a configuration with y(t) = hy a 
generic gauge transformation (I17p . Indeed, the energy fmictional is positive definite when y(t) = and 
negative definite when x{t) = 0. 
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Stokes line, the asymptotics of the solution to the differential equation changes its 
nature. 

• Another approach is to solve the system (fT8l) . Were the constraint G\E' = not 
imposed, the spectrum would be Emn = n — m with the eigenfunctions \nm) = 
\n')x\'>Ti)y. It is infinitely degenerate at each level. The constraint G = picks up 
only one representative of the set of eigenstates of H with a given energy. For 
example, the zero energy state annihilated by G is 

^0 = f:i2A:,2A:)(-l)'=^^|^. (20) 

At large k, the coefficient is proportional to l/v^, i-e. the normalization integral 
for (120|) diverges logarithmically. 

Similarly, only one eigenstate is left at each energy level. The full spectrum is 
discrete. 

En = n (21) 

with positive, negative, or zero integer n. It is unbounded both from below and 
above. This notwithstanding, the spectral problem is well defined and the evolution 
operator 

]C{x',x) = ^vl/:(a;')vl/„(a;)e*"* (22) 

n 

is unitary!^. 

Comparing the results we obtained under two quantization procedures, one can make 
two observations. First, the spectrum is shifted by 1/2. The ambiguity whether En = n 
or En = n + 1/2 has the same nature as the well-known ordering ambiguity — there are 
many different quantum problems having the same classical limit. The second observation 
is that, on top of the ordering ambiguity, there is in this case also another ambiguity 
associated with gauge choice. With any gauge choice, half of the spectrum involving 
either the states with negative or with positive energies is lost. 

A lesson that can be drawn from this simple toy model is that, for gauge systems, 
fixing the gauge classically and quantizing afterward may be dangerous. Certain essential 
features of the spectral problem (ITS]) may be lost. 

3 The potential —{ix)'^^^^. Classical dynamics. 

Having being equipped with necessary tools, we may proceed now with the analysis of 
the Hamiltonians of interest written in Eqs. p|7|) . We will concentrate mainly on the 
Hamiltonian ([7j) without the oscillator term in the potential. 

^See Ref. [23] for detailed discussion of this and related issues. 



9 



Let first n = 1. Consider the complex Hamiltonian 



n = — + iz^ (23) 



with z = X + iy,^ = p — iq. Its real and imaginary parts are 



G = -pq + x^ - 3xy^ . (24) 



y 
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Figure 3: Family of trajectories with E = 64. 

Consider the dynamics of the system described by the Hamiltonian H and the con- 
straint G. It can be treated as a gauge system. The equations of motion follow from the 
Hamiltonian H + AG, where A is the Lagrange multiplier. They have the form 

p = 6xy + 3\{y^ — x'^), x = p — Xq, q = 3{x'^ — y'^) + QXxy, y = —Q — Xp, 

G = -pq + x^ - 3xy^ = . (25) 

The Lagrangian (HM is invariant up to a total derivative with respect to gauge trans- 
formations f[T7l) with time- dependent parameter a{t). To find the classical solutions, we 
need first to fix the gauge. A convenient partial gauge fixing corresponds to the condition 
X(t) = 0, in which case the equations are reduced to 

p = 6xy, X = p, q = 3{x'^ — y^), y = —q, G = 0. (26) 

The solutions to (126|) belong to two classes: 1. Runaway trajectories, which reach 
infinity at finite time. These are, for example, the trajectories with initial conditions 
x(0) = x(0) = 0. They run away in the positive y directions. 2. Besides, there are 
families of closed orbits related to each other by gauge transformations (JT7I) with constant 
a. For positive energies, these families, depicted in Fig. 3, were found in Ref. [19]. This 
family has one distinguished member (one can call it a stem trajectory): the trajectory 



10 



which connects the turning points (the points where the monenta p, q vanish) with the 
coordinates. 



y* 



^1/3 



(27) 



(there is also the turning point x = 0,y = E^^'^, but the trajectories starting run away 
rather than coming back). 

Note that the famihes of trajectories with negative energies also exist (see Fig. 4). 
They stem from the trajectories connecting the turning points 



-E) 



1/3 



/3 



(28) 
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Figure 4: Family of trajectories with E = —64. 

Let us calculate for future purposes the action on these trajectories. Using the fact 
that the action for all orbits belonging to one family is the same, one can write 

{pdx + qdy) = 2Re / irdz = 2Re / ^/2{E - iz^)dz , (29) 

J Zl J zi 

where 2:1.2 are the turning points. For the trajectories of positive energies, the integral 
can be easily done by deforming the contour such that it passes the origin □, 

S+ = 4Re [ J2{E - iz^)dz = 4Re{z,)V2E [ ds = v/e^^s/e ^(4/3) ^ ^^^^ 

Jo Jo r(ii/6) 

To calculate the action for negative energy orbits, one has to take into account the fact 
that the turning points are at the same time the branching points of the integrand in 
(!29l) . For positive energies, the corresponding cuts do not hinder the deformation of the 
contour, but, for negative energies, they do. The cuts should be drawn such that the 
original path does not cross them. The deformed contour also should avoid crossing the 
cuts. The corresponding structure of the cuts, the original and deformed contour are 
shown in Fig. 5. 



^This result (in somewhat different normalization) was obtained in [19]. 
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Figure 5: Analytic structure of 7rci(-2) for negative energies. Solid line - original contour, 
dashed line - deformed contour, thick solid lines - the cuts. 

It is clear from the figure that the deformed contour involves four pieces: (i) from 
the left turning point to the origin, (ii-iii) from the origin down the cut and up again, 
(iv) from the origin to the right turning point. A simple analysis shows that the con- 
tribution of the parts (i-iv) involves an extra factor sin(7r/6) = 1/2 compared to the 
contribution of the parts (ii-iii). All together, the integral for involves an extra factor 
[l-|-sin(7r/6)]/ cos(7r/6) = -\/3 compared to the integral (|30i) for 5+ with the same absolute 
value of energy. In other words. 




Figure 6: Turning points, cuts and stem trajectories for the potential —iz^. Positive 
energies. 
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Figure 7: The same for negative energies. The cuts are drawn not to interfere with the 
paths. 

Again, there are runaway trajectories taking a finite time to reach infinity in the 
positive y direction. Besides, there are four famihes of closed orbits: two famihes with 
positive energies and two famihes with negative energies. The structure of the turning 
points, associated cuts and the stem trajectories connecting the turning points is shown 
schematically in Figs. 6,7 for positive and negative energies, respectively. Let us find the 
classical action at these trajectories. For positive energies, 

TT r(6/5) 



= 4cos-v^ rv^r^t/sE^/i" = 2v^cos-4^^^E^/i° 
+ 10 io 10 r(17/10) 

For negative energies, 

St"" = 2V^(l + sm5^)jI|5^(-B)V.o. (35) 



^"The sign of the potential corresponds to the convention ([7]) and to the conventions of [9, 19]. These 
conventions are convenient to make the physics of the systems ([7|) with different n more similar. 
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The superscript "up" in Eq. (!35|) refers to the upper trajectory in Fig. 7 going between the 
points e^^™^^^ and e~*'^/^°. The result for S'^^'^ is obtained in the same way as the result 
(13T!) . with the factor 1 + sin(7r/6) being replaced by 1 + sin(37r/10). When deforming the 
contour for the upper trajectory, we find, in addition to the parts composing the deformed 
contour of the lower trajectory and giving the factor 1 + sin(37r/10), also two extra pieces 
with the contribution ~ sin(37r/10) + sin(7r/10). The origin of all these factors can be 
clearly seen, if deforming the contour and the cuts in the way shown in Fig. 8. All the 
pieces (of nonzero length) connect the branching points to the center of the pentagon 
z = 0. 

Q 




Figure 8: The deformed upper contour wiggling around the deformed cuts. 

By the same token, for the potential — (iz)^""*"^, there are 2n families of the trajectories: 
n families with positive energies and n families with negative energies. 

As we have seen, the classical dynamics of the system with the potential — (iz)^""*"^ 
is similar in many respects to the complex oscillator dynamics: a distinct feature of 
both systems are the families of closed orbits with positive and negative energies, the 
members of one family being interrelated by gauge transformations. There are also two 
important differences. First, the system — (iz)^"^^ involves besides closed orbits also 
singular runaway trajectories. Second, for the complex oscillator, the stem trajectories 
for the families of orbits could be conveniently obtained by fixing the gauge y = or 
X = 0. But for the system —{iz)'^"'^^, this is not true. To begin with, the stem trajectories 
displayed above are essentially complex. This observation is not yet sufficient, however, 
because it does not exclude a conceivable in principle possibility that the trajectories can 
be put onto the real (or imaginary) axis by a complicated gauge transformation (ITTl) with 
nontrivial a{t). 

Let us find out what happens if we do fix the gauge y = for the system fl24f25p . 
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From G = 0, we deduce q = x^/p and hence the Hamiltonian is reduced to 

2^ 



2 fi 

H' = ^-lr2- (36) 



The corresponding equations of motion 



3x^ 

x=p+ — , P = ^ (37) 



follow from fl25p with A = —x^/p^. We see now that the reduced Hamiltonian 
is neither positive nor negative definite and involves only runaway trajectories. Closed 
orbits have disappeared ! This is another manifestation of the fact discussed in the 
previous section that fixing the gauge at the classical level is not an innocent procedure 
and may lead to a loss of important dynamic features. For the complex oscillator with 
the gauge choice y = 0, half of the orbits (the orbits with negative energies) were lost. 
For the system — (z^;)^""*"^, all closed orbits are lost and we are left only with runaway 
solutions. 

Let us discuss the relationship of the Hamiltonian (!36|) to another Hamiltonian ob- 
tained from (l23l) by a non-unitary rotation technique in the spirit of [15]. Let us multiply 
the potential by a coupling constant g, ix^ — ^ igx^, and find an operator R such that 
the rotated Hamiltonian H = e^ijP' /2 + igx^)e~^ be manifestly real. Then R can be 
presented as an infinite series over the coupling constant, 



and [24] (see also sect. V of Ref. [25]) 

^ = ^ + ^ + 0(/). (38) 
2 ap^ 

We see that the H* and H have similar structure, but the coefficients differ. This does 
not represent a paradox because H, in contrast to H*, involves the whole infinite series 
in g. Anyway, all the terms in this series are nonlocal, and one cannot obtain from this, 
say, the spectrum of quantum Hamiltonian as a perturbative series in g. The nonunitary 
rotation techniques is better suited to the problems like (jl]), where all the terms in the 
perturbative series for H are local. 

Coming back to fixing the gauge with the condition y = 0, it does not work well also 
for the mixed system (jl]), however small g is. The extra piece in H* is still nonlocal and 
singular at the turning point of the unperturbed oscillator trajectory where momentum 
p vanishes. As a result, the trajectory does not turn there, but rather stumbles and runs 
away. 



4 Quantum dynamics. 

Let us discuss now quantum dynamics of the Hamiltonians fl23l) . f l32|) . Consider Eq. fl23l) 
first. In Sect. 2, we outlined two regular ways to quantize gauge systems: ('z^ resolving the 
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constraint (s) at the classical level and quantizing afterward, and (ii) solving the system 
of differential equations (ITS!) with proper boundary conditions. 

To resolve the constraints classically, one has to fix the gauge. Unfortunately, as we 
have just seen, it is difficult to find a clever way to do it in our case. A natural gauge 
fixing leads to the problem involving only runaway trajectories. This means trouble and, 
indeed, for the highly nonlocal and not positive definite Hamiltonian fl36l) . one cannot 
formulate a well-defined quantum problem with a unitary evolution operator. 

Another approach is to solve the system f|T8|) . This is a nontrivial numerical problem. 
Indeed, one- dimensional spectral problems can be easily solved with Mathematica, but in 
this case the problem is essentially two-dimensional, which is much trickier. What is even 
more important, the operators H and G in f |T8l) are not elliptic, as usual, but hyperbolic. 
It is not thus evident that a reasonable solution to this problem exists... We will discuss 
this question somewhat more in the last section, but, basically, we leave it for future 
studies. 

There is, however, a way to define a consistent spectral problem related to the Hamil- 
tonian fl23l) [9] . Forget for a moment all what was said above about complexification and 
consider the Schrodinger equation at the real axis. 



2dx^ 



9 . 3 



^ = . (39) 



with the condition that the wave function falls down at a; = ±cxd. It is convenient to pose 
the problem not on the whole line (— oo, oo), but on the half-line (0, cxo). One can do it by 
exploiting the PT-symmetry of the potential (the property V{~x) = V*{x)). It dictates 
that for any solution \E'(a;) of Eq. flS^ . the function ^*(— x) is also the solution with the 
same eigenvalue. The functions 

^^(x) = ^(x) ±^*(-a;) (40) 

with the symmetry properties "^-^{—x) = \E'!j_(x) and \E'_(— x) = — \E'l(x) also satisfy this 
equation. We are hence allowed to consider the equations for PT-even function \E'+(x) 
and PT-odd function '^_{x) separately. In this case (in contrast, e.g. to the standard 
oscillator problem), the equation for \Ef_(x) does not give anything new. Indeed, one can 
make a PT-odd function out of a PT-even one by simply multiplying the latter by i. A 
generic solution to fl39|) is obtained by multiplying a PT-even solution by an arbitrary 
complex factor. 

The condition \E'(— x) = '$*{x) means that ^'(O) is real while \E''(0) is imaginary. 
By turning computer on, everybody can be convinced that the equation fl5^ with the 
boundary conditions 

Re f 1^1=0, ^{oc) = (41) 



^(0) 

has solutions at real positive discrete values of E. The remarkable fact is that these values 
are very close to semiclassical energies associated with the family of the closed orbits in 
Fig. 3 obtained from the quantization condition 

SiEk) = 7T{2k + l) , (42) 
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with the function S{E) being given by Eq. (!30|) . When k — oo, the spectral values 
extracted from Eqs. (13911411) and the semiclassical values extracted from Eqs. fl30f42l) rapidly 
converge. The exact and semiclassical values for Ek for first few levels [9] are shown in 
Table 1. 



k 





1 


2 


3 


exact 


0.763 


2.711 


4.989 


7.465 


scmicl 


0.722 


2.698 


4.980 


7.458 



Table 1. Exact and semiclassical spectra for the potential ix^. 

Once the solution is obtained, one need not to stay on the real axis. Actually, the 
solution can be continued analytically to complex values of the argument z in the regions 



Arg(2;) 



TT 

To 



< 



TT 

5 



Kig{z) 



llvr 



10 



< 



TT 

5 



(43) 



In other words, the spectral problem 

1 d 



2<l>2 



3^3 



E^ 



^ fd^/ds 
Re I 



0, ^(oo) = 0. 



s=0 



with $ = e'", still has a solution when a lies within the interval 
values are exactly the same as for the problem fl39ll4ip . When 

- 77r/10 <a< -37r/10 , 



(44) 

and the spectral 
(45) 



the spectrum is continuous: any positive or negative energy is acceptable. This is espe- 
cially clearly seen for a = — 7r/2 (meaning $ = —i). The problem is then reduced 
to 



1 d 



E^ 



^(oo) 



2ds^ 

0, Im [^'(0)] = 



(46) 



The real part of \E''(0) is not fixed, however, and tuning this parameter, one can obtain 
the solution dying at infinity at any energy, [jj A numerical analysis shows that it is true 
in the whole interval fHSj) . 

On the other hand, for tt/IO < a < 97r/10 the problem (jH]) has no solution whatsoever: 
the spectrum is empty. This is all illustrated in Fig. 9. 



""^^By modifying the spectral problem by, for example, imposing the conditions ^'(0) — '^{oo) = 
instead of ([1T|) . one can force the spectrum to be discrete and negative definite. But the condition 
^"(0) = is artificial and has no physical motivation. In particular, the discrete negative definite spectrum 
thus obtained has nothing to do with the semiclassical spectrum (|47)) . 

If lifting the requirement that the wave function dies away at infinity, the spectrum would again 
become continuous. 
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Figure 9: Spectral problem (jS]) in the complex z plane. 

The system of the lines separating the sectors in Fig. 9 form together with the positive 
imaginary axis the system of the asymptotes of the Stokes lines of the Schrodinger equation 
with the potential iz^ . Generically, for a polynomial potential of order n, such system 
involves n + 2 lines forming equal angles 2-7?/ (n + 2) [22]. 

The spectral problem (jUl) corresponds to the family of the classical orbits in Fig. 3 
with positive energies. As we have seen (in Fig. 4), there are also orbits with negative 
energies. Using the result (I3T!) . it is not difficult to find the corresponding semiclassical 
energies, 



Eh 



(2fc + i)r(ii/6) 
r(4/3) 



6/5 



TT \3/5 



1^ 



(47) 



One may suggest that a spectral problem should exist for which Eg. lHTl) would represent 
a semiclassical approximation. However, no such problem is known. |^ At least, it is not 
known in the standard form of boundary problem for some differential operator. One still 
can calculate the "exact spectrum" of such nonexisting (or very well hidden) problem by 
calculating corrections to the result (H7j) and representing -Eexact a series in semiclassical 
parameter ~ l/S'd. As this series is probably asymptotic, this method gives an intrinsic 
uncertainty in the spectrum ~ exp{— C^d}. However, the closeness of exact energies 
of positive energy states and their semiclassical approximations (see Table 1) and the 
calculations of higher order corrections in [19] suggests that this uncertainty is not large 
even for the "sky state" in Eq. fHTl) with k = and Sd = t^- It rapidly decreases with 
increase of k. 

Consider now the Hamiltonian (!32l) . Again, one can solve the Schrodinger equation 



^^And here is an important difference with the complex oscillator problem discussed in Sect. 2, where 
the spectral problem with the spectrum Ek = —k — 1/2 was perfectly well defined. 
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with the potential —ix^ at the real axis with boundary conditions \l'(±oo) = and find a 
discrete spectrum with real positive energies. As is seen from Table 2, these exact energies 
are very close to semiclassical values determined from the quantization condition 

^7 = 7r(2A; + l) , (48) 

where given in Eq. flMl) is evaluated for the upper trajectory in Fig. 6. 



k 





1 


2 


3 


exact 


0.710 


2.660 


5.458 


8.788 


semicl 


0.543 


2.608 


5.410 


8.750 



Table 2. Exact and semiclassical spectra for the potential —ix^. 

As we see, semiclassical approximation works somewhat worse in this case than for 
the potential ix^. But it works. 

We can now leave the real axis and solve the spectral problem 



d 



Re 



2$2 ds^ 
r d^/ds \ 



= 0, ^(cx)) = , 



(49) 



with $ 



e*". The solution with the same spectrum exists for 

137r 



a 



14 



For 



3n 

a -\ 

14 



TT 

^7- 



TT 
< -. 



or 



a 



7' 



or 



14 



177r 



a 



14 



TT 

< - 

- 7 



TT 

< 



7 



(50) 



(51) 



the solution still exists, but the spectrum is different. Its semiclassical approximation 
comes not from the quantization condition (l48ll . but rather from the quantization condi- 
tion 



S 



down 



7r(2A; + 1) 



(52) 



derived for the lower stem trajectory in Fig.6. The exact and semiclassical energy values 
for this case are given in Table 3. 



k 





1 


2 


3 


exact 


1.163 


5.234 


10.795 


17.428 


rpk 
semicl 


1.080 


5.186 


10.759 


17.400 



Table 3. Exact and semiclassical spectra for the potential —iz^ in the region (15T|) . 
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Figure 10: Spectral problem (H9ll in the complex z plane. 

Finally, for — 97r/(14) < a < — 57r/(14), the spectrum is continuous while, for 37r/(14) < 
a < ll7r/(14), the spectrum is empty. The corresponding regions in the complex z plane 
are displayed in Fig. 10. 

For the Hamiltonian (jlj) , there are n different nontrivial spectral problems with discrete 
spectrum defined in the sectors 



[In — l)7r 27rm 



2(2n + 3) 2n + 3 



< or mirror image , (53) 

- 2r2 + 3 ^ ' ^ ' 



m = They correspond to n different families of classical orbits with positive 

energies for the potential — flz)^""*"^. The problem studied in details in Refs. [9,19] was 
defined in the sector m = 0. lj We concentrate in this paper on imaginary potentials like 
in Eq.Q. But one can equally well [9,19] consider the potentials 

V{z) = -{izf'' (54) 

For the quartic potential ~ —-2^, there are four turning points, two sets of symmetric 
classical positive energy orbits and a corresponding spectral problem defined in the sector 
|a — 7r/6| < 7r/6 and its mirror images. For the potential ~ z^ we have, besides the 
standard spectral problem on the real axis, also a nontrivial problem in the sector \a — 
7r/4| < 7r/8, etc. For generic n, the potential fl54l) admits n/2 different spectral problems 
when n is even and (n + l)/2 different spectral problems when n is odd. 



"'^^There are also asymmetric spectral problems. One can, for example, go from infinity to zero along 
the line a = 177r/14 and to infinity from zero along the line a — tt/IA. But such problems have complex 
eigenvalues [20] and we are not considering them. 

In Ref. [19], the problem with m = 1 was also considered. It was represented as the problem with 
the potential V{x) ~ x^{ixY. But the results for the spectrum were given there only for negative e. 
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The presence of several different quantum problems associated with a given classical 
potential seems to be natural in view of our analysis for the complex oscillator, where two 
different spectral problems exist. However, it might appear surprising in the framework 
of Mostafazadeh's approach where the crypto-Hermitian Hamiltonian is obtained by a 
nonunitary rotation out of Hermitian H representing a quite definite series in g and hence 
the spectrum of H and of H represents a quite definite series in g. For example, the 
ground state energy of the system 



H 



igx 



(55) 



IS 



En 



2 32 ' 



(56) 



The resolution of this paradox is the following. Seemingly, only one of the spectral 
problems fH9l) associated with the Hamiltonian fl55|) . the problem defined in the sector 
including the real axis, can be safely treated in the framework of Mostafazadeh's approach. 
The ground state energy is plotted in Fig. 11 as a function of g. Indeed, the spectrum 
tends to the oscillator spectrum when g ^ It is not seen on the plot, but for very 
small g , starting from g — .02, the numerical values of the energies agree with the 

perturbative evaluation (l56l) . 



Energy 




0.2 0.4 0.6 



Figure 11: Ground state energy of (!55|) as a function of g. Upper region. 

The solution for another spectral problem at the vicinity of the rays a = —Sir/IA and 
a = 177r/14 behaves in a different and rather unexpectedly different way. For very small 
g, the spectrum is transformed, indeed, to the oscillator spectrum, but this transformation 
occurs in a very nontrivial manner. When g goes down, the energies of all the states go 
down in such a way that the energy of the ground state gets closer and closer to the 
energy of the first excited state. At some critical value of the coupling g^ ^ .03717, their 
energies coincide, 

Eo{g,) = E,ig,) ^ .484. 
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At still lower values of g, the energies should become complex. On the other hand, the 
second excited state goes down and down with decreasing of g and approaches the ground 
state oscillator energy without adventures, E2{g — > 0) ^ 1/2 (see Fig. 12). 

Energy 




log g 



Figure 12: Three first levels of (l55ll as a function of g. Lower region. 

To be more precise, there are no adventures in a sense that there is no phase transition 
and the state exists at any g and has real energy. But the asymptotics is reached only at 
rather small couplings. The energy of the second excited state finds itself at the vicinity of 
E = 1/2 only at (7 ~ .01. Now, i?2(-01) ~ .46 and does not coincide with the perturbative 
expansion (!56|) . It is not excluded that at still smaller values of coupling, g ^ .001, the 
perturbative asymptotics (!56l) finally shows up. To see whether it is true or not, a more 
careful numerical study is required. 

The third and the fourth excitations of the Hamiltonian ( l55l) coalesce and their energies 
cease to be real at g^^ ~ .007 (the energy is E^^^ ^ 1.37 at this point), while the fifth 
excitation approaches the first oscillator excitation E = 3/2 at very small values of g. 
One can suggest that this pattern holds also further up: the 6-th and the 7-th excitations 
of the mixed Hamiltonian coalesce and their energy becomes complex at some very small 
(7***, while the 8-th excitation approaches the second oscillator excitation E = 5/2, and 
so on. We thus observe an infinite sequence of "phase transitions" in the coupling, o 

This analysis shows that the Hamiltonian (155|) is crypto-Hermitian for all couplings 
in the upper sectors in Fig. 10, but, in the lower sectors, it is true only for not too 



This kind of transition when a pair of real eigenstates of a boundary problem coalesce and become 
complex is a known phenomenon [7]. Its essence is clearly seen in a trivial example. The matrix 

1 1 



a 1 

has a pair of close real eigenvalues at small positive a and a pair of complex conjugated eigenvalues for 
a < 0. When a — (the exceptional point), the matrix represents a nondiagonalizable Jordan block. 

An infinite set of such transitions in the parameter e for the problem V{x) = x'^{ixy was observed 
in [9]. We observed a similar phenomenon in a completely different physical context: it happens that 
some domain wall solutions in supcrsymmetric gauge theories disappear when mass of the matter fields 
exceeds certain critical values [26]. 
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small g > (7*. When g < g^,^ a. pair of complex conjugate eigenvalues should appear. For 
g < g-y^,, there are two such pairs, etc. It would be very interesting to see these complex 
eigenvalues explicitly. Unfortunately, it is not so easy to do it with our methods - the 
spectral problems of the type (HH HH]) make sense only for real energies — the boundary 
condition 



was derived under the assumption that "^{z) and \1'*(— 2;) satisfy the same Schrodinger 
equation, which is only true when E is real. A special study of this issue is required. 

5 Discussion and Outlook. 

Crypto- Hermitian systems have many common features with the systems involving higher 
derivatives. In both cases, Hermiticity of the Hamiltonian and unitarity of the evolution 
operator seem to be lost, but, if treating the problem properly, it is often restored. There 
exist also a more concrete relationship between two kind of systems. We have seen that the 
real part H of the complexified Hamiltonian [see Eqs.([9]), fl2^ ] is never positive definite 
and may give rise to ghosts. The same is true for higher-derivative theories. Actually, 
the canonical Hamiltonians of the latter have a rather similar form with not positive 
definite kinetic term [23] . The resemblance between the supersymmetric system analyzed 
in Ref. [23] and the problem considered here is even more striking. A system of the type 
fl24l) involves besides H the integral of motion G, and we are interested with the sector 
G = 0. The system studied in Ref. [23] (the bosonic part of its Hamiltonian is 



where (p, x) and (P, D) are two pairs of canonic variables and superpotential V{x) is an 
arbitrary function) also possesses an extra integral of motion N = P^/2 — V{x). In the 
sector with a particular value of N (including = 0) , the spectrum is discrete involving 
positive and negative energies. 

The latter is also true for the spectrum fl2T|) of complexified oscillator when the con- 
straint G = is imposed on the quantum states as in Eq. ffTSj) . The Dirac quantum 
problem flTSl) is more naturally posed than other quantum problems associated with the 
classical system in hand. This problem is easily solved in the oscillator case, but, for the 
potential iz^, this is a difficult immerical problem, and we leave it for further studies. One 
can speculate that its spectrum involves positive and negative energies, as the spectrum 
of the complexified oscillator and the spectrum of the Hamiltonian (!57|) do. However, 
it is an open question at present whether the problem (fT8l) makes sense for potentials 
more complicated than z^. As we have seen in [23], the Hermiticity of the Hamiltonian 
(1571) and the unitarity of the corresponding evolution operator are corollaries of the fact 
that classical trajectories of this system are benign enough: there are no collapsing or 
runaway trajectories where a singularity is reached at finite time. On the other hand, 
for the systems (jl]), runaway classical trajectories exist. For sure, not all the trajectories 
associated with the systems (jl]) are runaway trajectories. There are also closed orbits. 




H = pP- DV'{x) , 



(57) 
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and a hope that the problem (fT8|) is well defined is associated with their existence. The 
presence of runaway trajectories may spoil the brew, however. 

Runaway trajectories definitely spoil the brew for the quantum problems obtained by 
resolving the gauge constraint G = at the classical level. This procedure gives benign 
sensible Hamiltonians for the complexified oscillator. However, the Hamiltonian fl36p thus 
obtained is not Hermitian and unitarity is lost too. 

There are, however, Hermitian and unitary quantum problems associated with the 
Hamiltonians (jl]) and ([7]). For one of such problems corresponding to the potential — 
igx^ in the sectors below the real axis, we discovered a rather interesting and nontrivial 
phenomenon: when the coupling constant g is decreased, certain quantum states coalesce 
and disappear from the physical (real energy) spectrum. The number of such phase 
transitions is infinite, which reminds an infinite number of phase transitions in e for 
the potential x'^{ixy observed in [9]. Another phenomenon that comes to mind in this 
respect is the marginal stability curves in N=2 SYM theory and other supersymmetric 
systems [27]. When crossing these curves, quantum states may appear and disappear. 
However, the mechanism for this is quite different there. 

Let us make somewhat unusual conclusion listing again not the results obtained (that 
was done above), but rather the points which are not yet clear. 

1. It is not clear whether the spectral problem f|T8l) is well posed for the potential (jl]) 
and, if yes, what is its spectrum. Is the evolution operator unitary ? 

2. It is not clear whether one can formulate the spectral problems with discrete spec- 
trum in the dashed region in Figs. 9,10 by resolving the gauge constraint at the 
classical level with a clever gauge choice. 

3. It is not clear why, in contrast to the complex oscillator case, we have not found 
for the potential (jlj) a spectral problem involving only negative energy states (the 
cryptoghosts !) and related to the sets of classical orbits with negative energies. 
Can such problem be formulated ? 

The final remark is that crypto-Hermitian systems may prove to be something more 
than a formal mathematical exercise. They can bear relevance for physics. Our own inter- 
est to these problems stems mainly from their relationship to higher- derivative systems. 
And we believe (the arguments were presented in Ref. [28]) that the undiscovered yet 
fundamental Theory of Everything is a higher- derivative field theory {not string theory) 
living in higher-dimensional space-time. 

I am indebted to P. Dorey, L. Mezincescu, and D. Robert for useful discussions. 
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